Strong charge-transfer excitonic effects and Bose-Einstein exciton-condensate in 

graphane 
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Using first principles many-body theory methods (GW+BSE) we demonstrate that optical prop- 
erties of graphane are dominated by localized charge-transfer excitations governed by enhanced 
electron correlations in a two-dimensional dielectric medium. Strong electron-hole interaction leads 
to the appearance of small radius bound excitons with spatially separated electron and hole, which 
are localized out-of-plane and in-plane, respectively. The presence of such bound excitons opens the 
path on excitonic Bose-Einstein condensate in graphane that can be observed experimentally. 

PACS numbers: 73.22.-f, 78.67.-n, 71.35.Cc, 71.35.Lk 



Despite the short-live of grapheneQ and its derivatives 
the understanding and control of their properties rapidly 
approach a maturity [3, Chemical modifications by 
oxidation^, functionalization and dopingQ have en- 
hanced foreseen applications in nanotechnology. Recent 
synthesis of the fully hydrogenatcd-graphcnc (named 
" graphane" ) @ , which has been predicted to be a wide 
band-gap insulator with gap of 5.4 eV0,[il, 

adds to the 

portfolio of carbon-based structures for nanodevice ap- 
plications 0, 0]. The stability of the new 2D material 
has been analysed and two geometries (chair and boat) 
have been proposed^- However, there are still few open 
questions that need to be addressed: i) how much H is 
really incorporated in the samples? ii) does graphane in- 
herit negative electron affinity of hydrogenated diamond 
samples used for electron-emitters? Hi) what is the role 
of electron correlations in the band-structure and screen- 
ing; are electron-hole effects as important as for carbon 
nanotubes[§- 11| or small as in diamond 12 1? The present 
letter address those issues, in particular the last three, 
by means of first-principles calculations based on many- 
body Green function theory. Our results indicate the pos- 
sibility of having an excitonic Bose-Einstein condensate 
(BEC) by continuous pumping of excitons by light. Sim- 
ilar to the recent induced photoluminescence in graphene 
oxygenation [isj . the strong spatial localization of the ex- 
citons in graphane points to a highly efficient defect- 
induced luminescence in chemically (or doped) modified 
graphane. In our ground state calculations of the most 
stable chair conformation of graphane (Fig|T](b)) we used 
a pseudopotential plane- wave ap p roach |14lll5| within the 
LDA approximation to DFT[16|. While excited states 
have been treated using the state-of-the-art many body 
approach [H, [l3-[l^ . 



Before addressing the excitonic effects in the opti- 
cal properties of graphane including the formation of 



BEC, we emphasize that although the dynamical stabil- 
ity of the ideal graphane in the chair conformation (FiglT] 
(b)) has been recently shown [20|, the calculated phonon 
spectra agrees only partially with the measured Raman 
dataQ. In particular the measured D peak agrees well 
with the calculated Raman active Eg phonon at 1340 
cm^^ while the D' at 1620 cm~^ most likely originates 
from impurities as H vacancies. This means that there is 
a large defect density in graphane making it a hole-doped 
semiconductor with H vacancies acting as center for in- 
duced luminescence (see below). However at low dop- 
ing H vacancies do not affect significantly the electronic 
structure, so that ideal graphane allows to obtain a real- 
istic description of the electronic and optical properties of 
this system. On the other hand, the stability of graphane 
as hole-doped semiconductor, suggests the possibility to 
achieve a metallic phase under strong hole-doping by dc- 
hydrogcnation. Under this condition wc have observed 
a strong electron-phonon rcnormalization of the phonon 
frequencies around the BZ center (Kohn-anomalv) [isj . 
Similar to the case of doped diamond [2l|, this is a man- 
ifestation of a strong electron-phonon coupling pointing 
to the possibility of a high T^ superconducting phase. 

Electronic properties: the band structure of graphane 
(FiglT] (c)) is mainly dictated by the sp'^ hybridization 
of carbon orbitals (see Fig. [T] (b)). This causes a 
band gap opening (3.4 eV at the DFT-Kohn-Sham level) 
with respect to graphene. The classification of states 
is especially simple at the F-point. The top/bottom 
valence(Eg)/(Aig) corresponds to the C-C bonding a 
states, while two intermediate occupied bands corre- 
sponds to C-H bonding a states (Aig). The bottom of the 
conduction band is the anti-bonding A2M C-H a state and 
the anti-bonding E„ C-C cr state appears at higher energy 
(4.38 eV). The other bands correspond to free-electron- 
like states. Electron self-energy effects computed at the 
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FIG. 1: First Brillouin zone (a) and unit cell and bases of 
graphane in the chair conformation (b). Blue and white balls 
represent carbon and hydrogen atoms respectively, (c) Band 
structure in LDA (full line) and GW approximation (red cir- 
cles) and projected density of states of graphane. The zero 
indicates the vacuum level position. 



GTy-levellla llTl. Il9l| modify strongly the band structure 
(red dots on FiglT](c)): the fundamental energy gap at F 
increases from 3.4 eV to 5.4 eV, while the gaps at the high 
symmetry points M and K become 14.2 eV and 15.9 eV, 
instead of LDA values of 10.6 eV and 11.8 eV, in agree- 
ment vifith recent calculations Q . Moreover vife found that 
the vacuum level is at about 4.9 eV from the top of the 
LDA valence band (see FiglT](c)) so that graphane has a 
positive electron affinity both in LDA (1.44 eV) and GW 
(0.66 eV). 

Optical properties of graphane: In Fig[2] we plot the 
imaginary part of the macroscopic dielectric function 
eM(^) for a light propagating along the graphane plane 
in the ai-direction. To reveal the physical origin of differ- 
ent features in the optical spectrum we compare ImeM('^) 
calculated (i) without taking into account both the inter- 
electron and the electron-hole correlations, LDA-RPA, 
that is the random phase approximation (RPA) on top 
of the bare LDA band structure, (n) RPA using the 
GW quasi-particle spectra (GW-RPA) which neglects 
electron-hole correlations, and, finally, {Hi) from the full 
solution of the Bethe Salpeter equation (BSE) which ac- 
counts the excitonic effects (isl. [Tsl [l9| . 

The LDA-RPA absorption spectrum (green line in 
Figl2|) do not show any significant feature near to the 
LDA band gap (3.4 eV) because the corresponding transi- 
tions, being dipole-allowed, have small oscillator strength 
as the overlap between the top-valence (localized on the 
C-C bond,) and bottom-conduction (localized on the C- 
H bond) states is smah. Pronounced features in eM{i^) 
are present at higher energies. The peak about 8.5 eV 
is related to vertical transitions from the valence band 
to the conduction band at large wave vectors, while the 
peak at 10.2 eV corresponds to transitions near to the 
M point. The large structure at 11.2 eV arise from a 
Van Hove singularity near M corresponding to transi- 
tions from the valence band to out-of-plane states with 
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FIG. 2: Imaginary part of the macroscopic dielectric function 
for light polarized along the graphane plane. The dashed and 
dotted lines in the inset indicate the position of the GW gap 
and vacuum level respectively. For light propagating along 
the z axes the spectrum (not shown) is flat up to 6 eV. 



Kohn-Sham energies between 3.37 eV and 4.14 eV. The 
main effect of the quasi-particle corrections is a strong 
global shift of the absorption spectrum to higher energies 
(GW-RPA curve in Figl^]). Since the optical transition 
energy increases the amplitude of eM(w) is reduced to 
satisfy the f-sum rule. However, we show that electron- 
hole correlations dramatically modify the shape of eA/('^) 
(see BSE curve in FigH]). The "bulk" of the absorption 
spectrum nearly comes back to the original LDA-RPA 
position, but with a significant redistribution of the os- 
cillator strengths from the higher to the lower energies, 
and an appearance of a number of pronounced excitonic 
resonances that are related to electron-hole pairs with 
electrons in the delocalized out-of-plane states [l5|. 

In this work we concentrate on the most promi- 
nent physical effect of the electron-hole interactions in 
graphane, the appearance of bound excitons below the 
GW gap (see inset in Figl2|). These excitations, which 
are completely missing in RPA, are responsible for the 
UV absorption of graphane. The absorption spectrum 
for light propagating along x shows two bound excitons 
[B and C) with large binding energies of 1.6 eV and 
0.3 eV, respectively. The solution of the BS equation 
also reveals the existence of a dark exciton (A), 2 meV 
below the first optically active i3-exciton. The exciton 
C is related to transitions from the highest valence band 
to the first out-of-plane band and its wave function is 
delocalized out-of-plane[l^. In contrast, the spatial ex- 
tension of the strongly bound excitons A and B is very 
small, rex ~ 5.0 A, as it can be seen in FiglH] (b) and 
(d). These excitons are formed from the states of the 
double degenerate (at F-point) Eg valence band and the 
Pl2u conduction band (the hole states in these excitons 



3 



(b) 





Exc. A (xz Plane) Exc. A (xy Plane) 

(c) m _______ 





Exc. B (xz Plane) 



Exc. B (xy Plane) 



FIG. 3: 3D-Shape of the low energy excitons wave function 
for fixed position of the hole (green circle). 



belong to the lowest, E 



and the highest, Eg^''. bands 
respectively). Their wave functions reflect the C-H anti- 
bonding character of the states with the electron mainly 
localized on top of the H-atoms. Since the hole is lo- 
calized on the C-C bonds, the creation of such excitons 
corresponds to a charge transfer from the middle of the 
carbon plane to the side planes on top of the hydrogen 
atoms (see Fig 13] (a) and (c)). 

It is instructive to compare absorption spectra of 
graphenejiij and graphane. Electron- hole interaction is 
important in both systems. However, in the former it 
leads only to a redistribution of the oscillator strength 
resulting in the appearance of a resonant exciton at 4.5 
eV, while in the latter it gives rise to a strongly bound 
exciton at 3.8 cV where the spectrum of graphcne is com- 
pletely flat. This can be used as an additional experimen- 
tal fingerprint of the presence of hydrogenated platelets 
embedded in graphene. 

Considering the light propagating in the y-direction 
we find that the role of the strongly bound excitons is in- 
verted: A exciton becomes optical active, while B exciton 
is dark. Since the two excitons are nearly degenerate and 
have identical oscillator strengths for the corresponding 
polarization of light the absorption spectra remain al- 
most unchanged. The reason for this is the symmetry 
of Bloch states involved in the formation of the exci- 
tons: the electronic states have A2« symmetry, while the 
hole states belong to the 2D irreducible representation Eg 
with the eigenfunctions V'^^ and "(Ab'- A group theoreti- 
cal analysis [l^ shows that the only nontrivial dipole ma- 
trix elements are: {ipA2^, l^^l^^"*) ^^"^ (V'A2„ Ij/IV'b which 
explains why the A and B excitons are only visible for the 



y- and cc-polarization, respectively. For the light propa- 
gating perpendicular to the graphcne plane both excitons 
are dipole forbidden. 

The presence of H-vacancies with the strong spatial 
localization of the lowest A-exciton makes graphane to 
luminescence with high yield upon UV excitations, simi- 
larly to layered h-BN where a high luminescence yield is 
observed from a strong localized exciton [23j|. 

A striking feature of the excitons A and B is their 
large binding energy, which at least by an order of mag- 
nitude exceeds that in typical wide gap semiconductors. 
This becomes especially surprising if we note that both 
the conduction band and the Eg valence bands are 
nearly perfect parabolas over a large fraction of BZ, and, 
therefore, the excitons, despite their small radius, should 
be reasonably well described within the effective mass 
approximation. Using the standard method of invariants 
[2J] we derived the k • p Hamiltonian at F for the valence 
bands 

Hik) = ai{kl + kl) + [3\a,{kl - k^ + 2a,k,ky], (1) 



2m, 



15[, from 



and for the conduction band i^(k;) 
which we find the value of the reduced electron-hole mass 
flex = 0.29mo, which is very typical for most known semi- 
conductors (here mo is the bare electron mass, and for 
simplicity we model the two valence bands by a single 
parabola with an average mass). A resolution of the 
apparent paradox is an unusual form of the effective 
electron-hole interaction in graphane. In 3D dielectric 
media the effective interaction is obtained by the replace- 
ment e^/r — > e^/er, where e is the static dielectric con- 
stant. This simple recipe does not work in graphane that 
is a 2D dielectric, where the very notion of the dielectric 
constant makes no sense. In fact, an external electric field 
E induces a polarization P which is boimd to the plane: 
P(r) = a2DE{r)S{z), where a2D is the internal polar- 
izability of the 2D dielectric. Using the above relation 
in the Poisson equation for a point charge we get the 
following effective interaction potential (its 2D Fourier 
component) 



K//(q) 



27re^ 



|q|(l + 27ra2D|q|) 



(2) 



which is very different from the trivial renormalization 
of charge in 3D systems. The only parameter entering 
Eql2]is the polarizability a2D of a single graphane layer, 
which we extract from our ab initio results. Since our 
calculation are performed for a periodic stack of layers 
with sufficiently large inter-layer distance L, the 2D po- 
larizability is related to the actually calculated 3D polar- 
izability as follows a2D = L^^, where e is the dielectric 
constant of our 3D multilayer system. We can also esti- 
mate the exciton binding energy E^x by the variational 
solution of the Schrodinger equation with the effective 
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potential in Eql2] using the simplest trial wave function, 
ip ~ e"''/", where a is the variational parameter 25|. The 
result El, = 2.0 eV is very elose to the ab-initio value 
of 1.6 eV, while the optimized value of the variational 
parameter a = 7.7a.u. almost perfectly matches the ac- 
tual excitonic radius. Hence we clearly see that unusual 
strong binding of excitons in graphane is a result of a 
weak and specifically nonlocal 2D dielectric screening. 

Because of the small radius and the large binding en- 
ergy of excitons it is tempting to consider graphane as a 
potential candidate for a realization of the Bose-Einstein 
condensation (BEC) of optically pumped excitons. In- 
deed, using the standard expression for the degeneracy 
temperature Tq of a 2D Bosc gasfiij: 
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(3) 



(here M^x = nic + « 1.3mo is the excitonic mass) 
we find that at room temperature Tg = BOOA' the quan- 
tum degeneracy of excitonic gas in graphane is reached 
for the density n « 5 lO^^cm^^, which is two orders 
of magnitude smaller than -^^^ = 5.1 10^'^cm^^ (the 
later value sets an order of magnitude of the Mott criti- 
cal density for the transition from the gas of excitons to 
a gas/liquid of unbound electrons and holes). Another 
argument in favor of excitonic BEC in graphane is a spa- 
tial separation of the charges in the excitonic state - the 
hole is localized in the middle C-layer, while the electron 
on top of the H-layers. Therefore graphane looks quite 
similar to the coupled quantum wells (CQW) structures, 
which are very popular for experiments on excitonic BEC 
27, 28 1 . In CQW 27, 2^ such a separation leads to a 
very long excitonic life time that reaches microseconds. 
To estimate the radiative life time r of the excitons in 
graphane we adopted a 2D version of the approach of 
Ref.[2^, which lead to the following relation 

1 Ane^VLcP 

where fi, d and A are the exciton frequency, dipole 
matrix element, and the unit cell area, respectively. The 
result we find, r = 15 ps, is not very encouraging - the 
value is too small to unconditionally speak about BEC. 
On the other hand, it is still larger that the lifetime of 
cavity polaritons where BEC has probably been observed 
[sol ] . Clearly a detailed analysis of the thermalization ki- 
netics is needed to make more definite predictions. De- 
fects will also play a role in pinning the exciton con- 
densate. Nonetheless, we believe that BEC of strongly 
bound, small radius excitons in graphane is an interest- 
ing possibility which deserves to be carefully studied both 
experimentally and theoretically. 
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NUMERICAL DETAILS 



Our calculations of the ground state properties are 
based on density functional theory [J Q (DFT) within lo- 
cal density approximation Q (LDA) implemented in the 
pseudo-potential plane wave frameworkQ. The chair 
conformation has been simulated using a two dimensional 
triangular lattice with a basis of two CH units where the 
ionic potential has been modeled using Bachelet, Ham- 
man and Schliiter pseudopotentialsj^. The distance be- 
tween two adjacent layers has been set to ISA to simulate 
an isolated Graphane sheet. We chose a 90 Ry kinetic en- 
ergy cutoff for the electronic wave function expansion and 
a 16 X 16 X 1 Monkhorst-Pack[^ mesh for Brillouin zone 
(BZ) integration. Then, the optimized unit cell has been 
obtained minimizing the total energy as a function of the 
lattice parameter. At each value of the lattice constant 
the atomic positions {i.e. the internal degree of freedom) 
were fully relaxed to minimize the Hellman-Feynman|7[ 
forces. In the optimized structure the two CH units relax 
up-word and down-word respect to the ideal position of 
the graphene layer resulting in the formation of sp"^ C-C 
and C-H bonds. In particular the bond length are 1.51 
A (for the C-C bonds) and 1.11 A (for the C-H bonds) 
in agreement with previous calculations [1, Q. 

Phonon frequencies were calculated using the plane 
waves method within Linear Response Theory [l^l using 
a 30 X 30 X 1 Monkhorst-Pack fc-point mesh to compute 
the dynamical matrix. This allows us to achieve a very 
good convergence on the phonon spectrum (about 2% 
in the phonon frequency). A 10 X 10 x 1 (/-mesh in the 
BZ was used to interpolate the force constants for the 
phonon dispersion calculation. 

Finally in calculating the QP energies we used the GW 
approximation of Hedinflll. [l2l |. and the Bethe-Salpeter 
equation (BSE) has been solved to cope with excitonic 
effects in the absorption spectra 1^ |13[ . In particular in 
our calculation of QP energies we used a 30 x 30 x 1 fc- 
point mesh and 100 bands for computing the self energy. 
We take into account dynamical screening effects in W 
through the generalized plasmon pole model. The same 
fc-point mesh with 40 bands has been used to evaluate 
the dielectric function in RPA. The electron-hole kernel 
in the BSE and the dielectric function has been evaluated 
on a 30 X 30 X 1 k-point mesh. Moreover we used 100 
bands for the calculation of W and 20 bands to evaluate 
ej\/(w) since we were interested to the low energy region 
of the spectrum. 
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FIG. 4: Phonon spectra of ideal (solid line) and 1% hole- 
doped (dasched line) graphane in the chair conformation 



PHONON SPECTRA FOR PURE AND DOPED 
GRAPHANE 



In Fig|3] (solid line) we show the calculated phonon 
dispersion curve for a perfect graphane in the chair con- 
figuration. The phonon spectra presents three well sep- 
arated types of vibrations: the low-frequency acoustic 
modes (not shown), the intermediate- frequency optical 
modes involving the stretching of the C-C bonds, and 
separated by a gap of 1500 cm~^ (not shown), we find 
optical modes corresponding to C-H bonds. The presence 
of this large gap proves that the vibrations of H atoms 
are almost completely decoupled from the other modes. 

Comparison of our calculated spectra with available 
experimental Raman data, shows that the measured D 
peak agrees well with the calculated Raman active Eg 
phonon at 1340 cm~^ while the D' peak observed at 1500 
cm~^ most likely is related to impurity effects due to H 
vacancies. 

Finally, in Fig|4] we compare the phonon spectra of 
ideal graphane (solid line) with that of 1% hole-doped 
graphane (dashed line) obtained in the rigid-band ap- 
proximation. As we can see hole-doping causes a strong 
renormalization of the phonon spectra near to the F 
point. In particular the two degenerate optical modes 
related to the C-C stretching with H atoms moving in- 
phase with C atoms soften from 1168 cm~^ to 542 cm~^ 
while the two degenerate optical modes related to C- 
C stretching with H atoms moving out-of-phase with C 
atoms soften from 1333 cm~^ to 1221 cm~^. The large 
Kohn anomaly at 542 cm~^ suggest a strong electron- 
phonon coupling in doped graphane related to the C-C 
covalent bonds. This makes doped graphane a possible 
candidate for high Tc superconductivity. 
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SYMMETRY PROPERTIES OF THE WAVE 
FUNCTIONS 

In this section we perform an analysis of symmetry 
properties of the wave functions corresponding to the va- 
lence and conduction bands at the F point based on the 
group theory. This allows to explain the behaviour of 
the absorption spectra at low energy for different po- 
larizations and the nature of the two nearly degenerate 
excitons. 

Graphane in the chair conformation belongs to the D^j^ 
point group which is characterized by 12 symmetry oper- 
ations (with inversion) and 6 irreducible representations: 
two 2D representations Eg and E„ and four ID represen- 
tations Aig, A2g, Ai„ and A2„ (see table |T|. 





E 


2C3 


3^ 


i 


2^6 


3crd 


Basis functions 


Ala 


1 


1 


1 


1 


1 


1 


2,2 2 
X +y ,z 




1 


1 


-1 


1 


1 


-1 






2 


-1 





2 


-1 





{x^ - y'^,xy) {xz,yz) 


Aiu 


1 


1 


1 


-1 


-1 


-1 




A2u 


1 


1 


-1 


-1 


-1 


1 


z 




2 


-1 





-2 


1 








TABLE I: Character table for Dsd point group. 

Since the valence bands arc two fold degenerate at the 
r point they belong to one of the 2D irreducible represen- 
tation Eg or Fju- However as can be inferred from FiglS] 
(a) and (c), the wave functions corresponding to the va- 
lence bands are symmetric respect to the inversion. This 
is sufficient to conclude that the valence bands indicated 
with e^^^ and with e^j belong to the Eg irreducible rep- 
resentation. Moreover we indicate with ^/j^^ and V'^^ the 
corresponding wave function, with the label 1(2) indicat- 
ing that the wave function is antisymmetric (symmetric) 
under reflection respect to the YZ plane. 

The conduction band is not degenerate and it is char- 
acterized by a wave function antisymmetric respect to 
the inversion (see Figl5](b)). On the other hand this lat- 
ter is antisymmetric under refection respect to the YZ 
plane. Thus the conduction band (eyi2„) belongs to the 
A2U (see tabl(il]) irreducible representation. 

Now we focus on the symmetry properties of a general 
polar vector (x,y,z) as for example the dipole moment. 
It can be easily shown that the third component z trans- 
forms in agreement with the A2U irreducible represen- 
tation, while the in-plane components (x, y) transform 
in agreement with the E„ irreducible representation. In 
particular y is the partner of x. Then, when the light is 
directed along the Z axis the possible states coupled with 
the valence bands by dipole moment are contained in the 
direct product ^ Eg — £'„. This means that transi- 
tions from the valence to the conduction bands arc forbid- 



(a) (b) 



(c) 




FIG. 5: Projection on the XY ((a) and (c)) and ZY (b) planes 
of the wave functions -0^^ (r) (a), 'ipA2u{r) (b) and "0^^ (r) 
(c). The red- blue lines refer to negative-positive values of the 
wave function. 



den for light directed along the Z axis and both excitons 
A and B are not visible. On the other hand, when the 
light is directed along X or F, the possible states coupled 
with the valence bands by dipole moment are contained 
in the direct product E^xEg = Aiu+A2u+Eu which is a 
4D reducible representation of the symmetry group in the 
bases defined by the four wave functions: X'ip^j^\ V^^e ^ 

xip^ and yipE ■ now possible, starting from this 

bases set, to generate a set of wave functions symmet- 
rically adapted with the three irreducible representation 
appearing in the direct product Eu x Eg. In particular 
we found that the only wave function which transform 
in agreement with the A2U irreducible representation is a 
linear combination of xip^^^ and yipE^'- This means that 
the only possible transitions involving the wave function 
V'A2„ are: (V-As., klV'^j) and (■(/'A2„ lyl^^j)- 

It is now clear why the most bounded exciton related 

(2) 

to electron-hole pairs ('!/'A2„ > '4'e2 ) optical active for 
light polarized along Y, while the highest energy exci- 
ton related to electron-hole pairs {ipA2u j ''PeI ) visible 
for light polarized along X. However, due to the small 
overlap between ipA2u and ipEg wave functions the two ex- 
citons present a very small strength compared with that 
of the higher energy excitons. As a results, it does not 
reflect with evident features on the behaviour of eM(w). 
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ANALYSIS OF THE EXCITONIC RESONANCES 
IN THE GRAPHANE ABSORPTION SPECTRA 
OF FIG. 2 

The prominent peak at 15.4 cV is completely de- 
pressed, while the absorption is strongly enhanced at 
10.2 eV and 11.0 eV. The first peak at 10.2 eV is re- 
lated to the correlated electron and hole states with large 
wave vectors from the Eg bands and out-of-plane bands 
with QP energies 0.9 eV < e^^ ^ 2.1 eV. The "binding 
energy" is about 0.5 eV. The second resonant exciton at 
11.0 eV has a lower binding energy (0.2 eV) and is formed 
by electron-hole pairs from C-H valence states and out- 
of-plane states with QP energy close to 1.20 eV. The 
series of excitonic resonances between 8.30 eV and 8.90 
eV are related to the electron-hole pairs corresponding to 
transitions from the small wave vector Eg states to the 
out-of-plane states with QP energies 0.13 eV< e^jf < 
2.02 eV. Finally, the lowest resonant excitonic peak at 
6.46 eV (peak D in the inset in Fig. 2) originates from the 
mixture of the transitions from the highest Eg valence 
band to the A211 conduction band and the first out-of- 
plane band. Because of the free-electron character of the 
electronic states involved in all these electron-hole res- 
onances the corresponding excitonic wave functions are 
completely delocalized in the out-of-plane regions (see 
Figl(b)). 



for each transformation g of the symmetry group. 
D^<'{g) being the 2D matrix representing g in the Eg 
irreducible representation. 

Since the matrix Eql5] represents iV^ = 4 equations for 
the elements Hij{IC) of H{IC), these constraints allow one 
to construct H{IC). 

To determine the correct form of the Hamiltonian, 
we first note that the tensor operator K. can be decom- 
posed into two irreducible tensor operators, one (/C"*^" = 
kl + ky) belonging to the Aig irreducible representation 
and one {JC^^ = (fc^ — ky,2kxky)) belonging to the Eg 
irreducible representation {2kj.ky being the partner of 
~ ^y)- Furthermore. H{1C) can be decomposed on 
a complete set of linearly independent 2x2 dimensional 
matrices that transform according to those irreducible 
representations which are contained in the direct prod- 



uct Eg X Eg 



Eg + Aig + A2g ■ Starting from the four 



linearly independent matrices / (identity matrix) , ax, 
and dz (Pauli matrices) we can construct a set of linearly 
independent matrices Xi symmetrically adapted with the 
irreducible representation T = Eg, Aig, A2g. This can be 
done in terms of the standard group theory using the 
relation: 



D^^igno''^ {g)Dl,{g) 



(6) 



(a) 





where x runs on the four matrices {i.e. the identity and 
the three Pauli matrices) and Dlj{g) is the matrix which 
represents g in the irreducible representation F. Eql6] al- 
lows to obtain from x ^ matrix symmetrically adapted 
with the first row of the representation F. The corre- 
sponding partners can be constructed from: 





Exc. C 



Exc. D 



FIG. 6: 3D-Shape of the excitonic wave function for the ex- 
citon C (a) and D (b). 



k p HAMILTONIAN 

In this section we derive the k-p Hamiltonian H{IC) for 
the valence bands around the F point [JC being a second 
order tensor operator which depends from the compo- 
nents of the wave vector k) . To determine the number of 
the independent parameters in H(IC) we use the method 
of the invariants 1^ , which consists to impose the condi- 
tion: 



Xi 



Y^b^^{g)xW^\g)Dl{g) 



(7) 



b^^{g)H{g-'K.)b^^\g)^H{lC) 



(5) 



with i ^ 1. In this case we found that / and ay trans- 
form in agreement with Aig and A2g respectively. While 
dz transforms in agreement with the first row of Eg and 
ax is its partner. To determine the matrices symmetri- 
cally adapted with Eg we chose {kx — ky,2kxky) as basis 
for the construction of D^!>{g). 

At this point, for each irreducible representation F, we 
can construct the corresponding invariant (7r) as a prod- 
uct of symmetrically adapted matrices xf a-nd irreducible 
tensor components A^f: = X^iXf^f- 

In particular we obtain two invariants: 'T^'^^ = /(^^ + 
ky) and T^" = az{kl — ky) + 2axkxky but we do not find 
any invariant transforming with the A2g representation. 
This happens because at the second order the tensor op- 
erator /C does not contain any term belonging to the A2g 
representation (see table [l|. On the other hand, in the 
presence of a magnetic field along the Z axes {Hz), one 
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other term gL^H^ appears in the Hamiltonian (L^ be- 
ing the Z component of the angular momentum) which 
results in a third tensor operator related to the commuta- 
tor [kx,ky] and belonging to the representation. As a 
consequence a third invariant T^^" oc ay [k^ , ky] appears 
in the Hamiltonian. 

Thus, each liner combination of the two invariants 
T"^^" and T'^" satisfies by construction Eq. [51 so that 
the correct form of the k ■ p Hamiltonian is: 



Hik)^aiikl + kl)+l3[a,ikl 



ky) -\- 2&xkxky] 



(8) 



with two independent parameters a and /3. Diagonal- 
izing Eq|8] we find that the band dispersion in k space 
is spherically symmetric with E(k) = (a ± -I- ky). 

On the other hand for the conduction band we found the 
trivial result: iJ(/C) = ^i^. 

Moreover, from the calculation of a (a = — 1.31/mo), (3 
(/3 = — 0.49/too) a-iid me (me = O.SStoq) by interpolating 
the ab-initio band structure (Fig.l) we found the mean 
value of the effective mass of the two nearly degenerate 
excitons close to ^ex = 0.29 toq (™o being the electron 
mass). 
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